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FOREWORD 


This report documents the research performed under contract 
NAS5-21860 between 7 September 1972 and 6 October 1973. Dr. M. 
Griggs was the Principal Investigator, and Dr. C. B. Ludwig and Dr. W. 
Malkmus were co-investigators. 

The author is indebted to Dr. R. S. Fraser, the contract monitor, 
for many stimulating discussions, and to Mr. G. Hall for making some of 
the ground-truth measurements at the Salton Sea. 
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DETERMINATION OF AEROSOL CONTENT 
IN THE ATMOSPHERE FROM 
ERTS-1 DATA 


M. Griggs 

Science Applications, Inc. 

SUMMARY 

Significant results, relating the radiance over water surfaces to 
the atmospheric aerosol content, have been obtained. The results indi- 
cate that the MSS channels 4, 5 and 6 centered at 0. 55, 0. 65 and 0. 75 ym 
have comparable sensitivity, and that the aerosol content can be deter- 
mined within + 10% with the assumed measurement errors of the MSS. 

The fourth channel, MSS 7, is not useful for aerosol determination due 
to the water radiance values for this channel generally being less than 
the instrument noise. The accuracy of the aerosol content measurement 
could be increased by using an instrument specifically designed for this 
purpose. 

This radiance -aerosol content relationship can possibly provide a 
basis for monitoring the atmospheric aerosol content on a global basis, 
allowing a base-line value of the global burden of aerosols to be established. 
This base-line could probably be established more rapidly from satellite 
measurements than from a network of ground-based observations. In addition, 
this technique could possibly provide a method for monitoring the particulate 
emissions of the SST's, by making observations in the vicinity of flight cor- 
ridors, such as over the North Atlantic. It may be possible to look at the 
ocean through the flight corridor and alongside it and measure the difference 
due to the SST's in the aerosol content. 
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Further studies of the radiance-aerosol content relationship 
are needed to establish the global applicability of the results, and 
to confirm that the effects of sun glitter are minimal as indicated in 
this program. These studies could utilize future ERTS-1 and ERTS-B 
data at the existing test sites. Cooperation with other agencies such 
as EPA and NOAA, who make Volz observations, would allow the studies 
to be expanded to national and global scales. ERTS-1 data obtained over 
the USA since August 1972 should also be analyzed in conjunction with 
aerosol data from the turbidity network of Volz photometers operated 
by EPA. 

The contrast -aerosol content investigation show useful linear 
relationships in MSS channels 4 and 5, allowing the aerosol content to 
be determined within + 10%. MSS 7 is not useful due to the low accuracy 
in the water radiance, and MSS 6 is found to be too insensitive. These 
results rely on several assumptions due to the lack of ground -truth data, 
but do serve to indicate which channels are most useful. 
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1 . 


INTRODUCTION 


The scientific community (e. g. SCEP^ and SMIC^) has become 
increasingly aware in recent years of the importance of atmospheric aerosols 
and their optical properties in possible climate modification. The aerosols 
in the atmosphere consist of man-made and natural particles, and it is the 
man-made contribution due to combustion added to the natural (dust, sea 
spray, forest fires and volcanic dust) background that is generally con- 
sidered to be important in determining climatic changes. However, the 
man-made contribution on a global scale is quite small; estimates range 
from a negligible amount^ to about 6% of the natural background^. (Of 
course, on a local scale, e. g. the Los Angeles basin, man-made particles 
can far exceed the natural concentration. ) It may well be that global changes 
in the natural background are more important than man-made particles. 

With the increased cultivation of land and activities of man in arid areas, 
the background level of aerosols is likely to increase. 

(5) 

McCormick and Ludwig' presented evidence of a worldwide buildup 
of atmospheric aerosols which could increase the earth albedo resulting in 
a cooling of the earth-atmosphere system. This effect would counteract 
the postulated increase of temperature in the lower atmosphere due to the 
"greenhouse effect" of the increased C0 2 emissions by human activities. 

In fact, there has been a decrease in the mean annual air temperature since 
about 1945 at mid latitudes, suggesting that the aerosol pollution effect is 
greater than that of the C0 2 increase. However, the effects of aerosols and 
C0 0 are more complex than suggested above, so that their effects on climate 
are not readily predicted. For instance, Robinsom ' points out that the 
earth may self-regulate its temperature by the variation of cloud amount: 
the higher temperatures, due to the C0 2 "greenhouse effect", lead to a 
higher water content in the lower atmosphere, which may increase the 
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cloud amount; this increases the albedo, thereby decreasing the temperature. 
Robinson concludes there is no justification for forecasting a final equilibrium 
temperature due to an increase in COg content, until atmospheric models are 
significantly improved to include the cloud cover as a variable. 

In addition to the uncertainties in the climatic effects of COg, the 
cooling effect of aerosols suggested by McCormick and Ludwig may not be 
correct. Charlson and Pilat^, Atwater^ and Mitchell^ have shown that 
since aerosols absorb and scatter, they may produce warming or cooling, 
depending on the ratio of absorption to scattering. 

Thus, it is clear that considerably more work on the complex problem 
of modeling the atmosphere and on the optical properties of aerosols is needed 
before the long term effects of man-made pollution can be predicted. Since 
these problems will not be solved in the near future, it is important to initiate 
global measurements of aerosols on a continuous basis to monitor any changes. 

We have shown in earlier theoretical studies^) that it should be pos- 
sible to make satellite observations of the aerosol optical thickness of the 
atmosphere from contrast measurements of ground features, and from ra- 
diance measurements. The launch of ERTS-1 offered the opportunity of 
investigating these theoretically derived relationships with actual satellite data. 
For convenience in comparing data at different wavelengths we shall refer 
to aerosol content as well as aerosol optical thickness. The aerosol content 
is defined in terms of the Elterman 1964 model vertical aerosol optical thick- 
ness; i. e. , the aerosol content is given by the ratio (measured aerosol optical 
thickness at wavelength X/model aerosol optical thickness at wavelength X). 

The aerosol optical thickness can be related to the mass loading of the 
aerosols using the relationship given by Griggs^^^. 

These relationships were investigated using MSS data at two test sites, 
with ground truth being obtained with ground-based Volz photometer measure- 
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ments of the aerosol content. The two test sites were at San Diego, where 
radiance measurements over the ocean were conducted, and at the Salton Sea/ 
desert region, where radiance and contrast measurements were made. The 
test sites are illustrated with ERTS imagery in Fig. 1-1 and 1-2. 
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Figure 1-1. San Diego Test Site 

ERTS MSS 6 Data 1-18-73 
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Figure 1-2. Salton Sea/Desert Test Site 
ERTS MSS 5 Data 12-30-72 
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2 . 


APPROACH 


The approach to the investigation has been an empirical one based 
on theoretical calculations for model atmospheres. To make the computations 
manageable, certain approximations about several parameters, such as the 
aerosol size distribution and the underlying surface reflectance, have to be 
made. Hence, in the real atmosphere, model conditions are never realized, 
so that deviations from the theoretical relationships are expected. Thus, an 
empirical investigation has been conducted using the theory to provide insight 
into the extremes of values which may be encountered. 

The two relationships studied have been one between radiance over 
water surfaces and the aerosol content, and another between the water/desert 
contrast and the aerosol content. The satellite radiance measurements were 
obtained from the ERTS-1 digital data, and the ground-truth measurements 
of the aerosol content were made with a Volz photometer at the time of 
selected ERTS overpasses. 

The theory and its limitations are discussed in the following sections. 

2. 1 Relationship of Radiance and Aerosol Content 

Calculations of the radiance backscattered from the earth-atmosphere 

(12 13) 

system, as seen from space have been published by Plass and Kattawar. y 
These calculations, using Monte Carlo techniques, consider multiple scattering 
of all orders, and take into account aerosol scattering and ozone absorption. 
Examination of the results of Plass and Kattawar shows that the outgoing 
radiance varies with aerosol content, and is most sensitive when the under- 
lying surface albedo is low. The ocean, which covers much of the earth, has 
a low albedo at high sun angles and provides a suitable underlying surface for 
aerosol measurements. The calculations also indicate that the longer wave- 
lengths are more sensitive to aerosol changes. At shorter wavelengths the 
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Rayleigh optical thickness is comparable to, or greater than, the aerosol 
optical thickness, so that changes in the aerosol content have less effect. 
The results obtained by Plass and Kattawar for 0. 7 fim and zero albedo 
are plotted in Fig. 2-1. 

In order to further investigate the effect of aerosols on the upward 

( 14 ) 

radiance over a calm water surface, Plass and Kattawar made some 
special calculations for us under our Contract NAS1-10466, for several 
aerosol vertical distributions. The calculations showed a result of great 
importance for satellite observations of the upward radiance: the upward 
radiance depends strongly on the total number of aerosols, but not on their 
vertical distributions. Thus measurements of the upward radiance can be 
directly related to the total vertical aerosol content and hence the global 
loading. 

The three wavelengths (0. 7 |im, 0. 9 ym and 1. 67 urn) considered 
show comparable sensitivity to aerosol changes. However, the relative 
normalized radiance is less at the longer wavelengths, and since the in- 
coming solar flux decreases at longer wavelengths the absolute radiance 
level decreases rapidly with increasing wavelength (the absolute radiance 
at 1. 67 ym is about 6% of that at 0. 7 (im). 

The results for 0. 7 jjm are plotted in Fig. 2-2 to show the relation- 
ship between the upward radiance, normalized to unit incident solar flux, 
and the aerosol content of the atmosphere for various sun angles. A simple 
linear relationship is shown to exist between radiance and the aerosol con- 
tent. These straight lines are based on only two values of aerosol content, 
but a linear relationship may be established by considering the four data 
points for zero albedo in Fig. 2-1. From these curves a knowledge of the 
sun angle and the absolute radiance at 0. 7 urn allows the aerosol content 
of the atmosphere to be determined. From Fig. 2-2 it is seen that a 1 

percent change in radiance is equivalent to about 1. 5 percent change in 

<! 

aerosol content. 
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The data in Fig. 2-2 may be used to compute the radiance in the 
MSS channel 6 , centered at 0. 75 4 m, by interpolation of the results of 
Plass and Kattawar at wavelengths 0. 4 4 m, 0. 7 4 m, 0. 9 4 m and 1. 67 4 m, 
and by assuming a rectangular spectral response 0. 1 4 m wide. The 
results, which neglect the oxygen absorption in this band, are shown in 
Fig. 2-3. 


2. 2 Relationship of Contrast and Aerosol Content 

(15) 

The theory of contrast reduction developed by Duntley is briefly 
discussed below. The inherent contrast of an object relative to the back- 
ground is defined as 


C 


o 



(2-1) 


where B and B ' are the radiances of the object and background, 
00 

respectively. The apparent contrast, as viewed from range, R, is 
defined as 


C 


R 




( 2 - 2 ) 


where B„ and B D ' are the corresponding radiances observed at R. 

R K 

Duntley shows that an object of radiance B q viewed through a scattering 
atmosphere with an optical depth, r, has an apparent radiance, B R , given by 
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Radiance (mw/cm ysr) 


0.3 


0.2 



ON 3N 

Aerosol Content (N = Elterman 1964 Model) 


Figure 2-3. Radiance vs. Aerosol Content for MSS 6 
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[l-e" T ] , 


(2-3) 



B (0) 
e _T + -S- 


in which the second term on the right results from scattering into the 

line of sight, where B (0) is the scattered radiance from unit thickness 

£1 

of the atmosphere near the object, and a is the atmospheric extinction 
coefficient near the object. 

Hence the contrast ratio C r /C q is found to be 



. B a<°> 

1 + a B' 
o o 


(e r - 1) 


(2-4) 


We fitted this equation with values of C r /C q calculated for the 
results of Plass and Kattawar^ 12 ^ and the values of B a (°)/ a 0 B 0 # obtained 
were strongly dependent upon r. The inapplicability of Eq. 2-4 to the 
case of observations made through the entire atmosphere is not unexpected, 
since the variation of r is due to the particulate component only while 
the Rayleigh component remains constant. 

A perfectly arbitrary generalization of Eq. 2-4 is made: 



1 


1 + g(r> A) 


(2-5) 


where A is the underlying surface albedo (assumed Lambertian). 

From values published by Plass and Kattawar for reflected radiance versus 
albedo, it was found that the reflected radiance is well represented by a 
linear function of A: 


R ref 1< T > A > - R refl< T ’ °> + A t R r C il (T - 1] ‘ R ren< T > °» • < 2 ' 6) 
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The apparent contrast is, consequently: 


^ _ refl' ’ ’ refl w ' 

C R= R en (T'7A-') 

(A-A')[R reI1 (r,l)-R refl {T,0)] 

= ^IrTTJTrA-Tftr ' ( 

where A and A ' are the albedos of the object and background respectively. 
Since 

C 0 - T < 2 - 8 > 

we have 



R refl (r ’ 0) 


(2-9) 


By identifying with Eq. 2-5 we find that the unspecified function g(r , A) 
may be separated into r - and A-dependent parts: 

g(r,A) = f(r)/A (2-10) 


where 


f(r) 


^refl 


^refl 


(r, O) 

- R 

refl 


( 770 ) 


( 2 - 11 ) 


Values f(r) using Eq. 2-11 have been calculated for the center wavelengths 
of three of the four MSS channels, based on interpolation of values of 
r( t , 0) and R(r, 1) calculated by Plass and Kattawar, ^ 12 » 13 » 14 ^ for the 
wavelengths 0.4 ^m, 0.7 ^m, 0.9 j^m and 1.67 jjm. The results, for 
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three sun angles, are shown in Figs. 2-4, 2-5, and 2-6. The fourth MSS 

channel between 0. 8 Jim and 1. 1 ym is not considered due to the strong 

absorption by water vapor, which is an atmospheric variable, in this band. 

The function f(r) is related to the experimentally observable quantities 

C , C D and A' by Eq. 2-5 and 2-10: 
o R 

1( T ) = A ' (C Q /C R - 1) . (2-12) 

Thus, a measurement of the contrast ratio and knowledge of albedo and 
sun angle yields a value of f(r) that determines t . 

The preceding analysis is based on calculations using a Lambertian 
underlying surface. Of course, in practice this type of surface does not 
exist, and reflectivity should be used instead of albedo in determining f(r) 
from Eq. 2-12. The reflectivities of natural surfaces vary with wave- 
length and sun angle and are not accurately known. A relatively simple 
case is the Salton Sea/desert contrast. However, the reflectivity of sand 
which has been measured for certain conditions by Coulson' ’ ; cannot 

be deduced for all the conditions required for this program. The radiance 

(18) 

reflected from water is more readily calculated, ' but varies with sur- 
face conditions and suspended matter. Thus, due to the shortcomings of 
the available reflectance data, and to the approximations in the theory, 
the empirical approach must be used. 

2. 3 Potential Problem Areas 

The discussions in Sections 2. 1 and 2. 2 are based on theoretical 
calculations which use a model atmosphere, model aerosol properties, 
and assume a smooth water surface or a Lambertian surface. Since these 
model conditions are never realized in practice, empirical relationships 
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f(r) 





between the satellite data and the aerosol content have been investigated. 
However, two potential problems must be considered, although as shown 
later in the discussion of the results, they do not appear significant. 

2. 3. 1 Sun Glitter 

If the ocean were perfectly smooth as assumed in the calculations 
an image of the sun would be seen at the specular reflection angle, and the 
only upwelling surface radiation observable at other look angles from space 
would be the diffuse sky radiation reflected from the ocean surface and the 
radiation scattered up from below the ocean surface. As the smooth ocean 
surface is increasingly disturbed, a glitter pattern becomes increasingly 
larger about the specular point. At sun zenith angles greater than about 
30° the glitter effect has been considered negligible (except for very rough 
seas) at the nadir point. However, measurements by Hovis (private com- 
munication from R. Fraser) suggest that this assumption is not correct, so 
that the ocean surface radiance at the nadir is not known accurately. 

This problem might be overcome by making observations at two 
wavelengths, assuming that the spectral variation of the surface radiance 
is known. The choice of wavelengths must be made carefully since the 
spectral distribution of the radiance does vary due to ocean properties 
such as chlorophyll content, suspended matter and depth. 

2. 3. 2 Surface Reflectance Gradients 

The calculations of Plass and Kattawar which are used in deriving 
the contrast attenuation relationship in Section 2. 2 assume an underlying 
surface of constant reflectivity extending to infinity. However, for the 
contrast measurements there are two adjacent surfaces of different re- 
flectivities, and radiation reflected from one surface is scattered into the 
atmosphere above the other surface so that the apparent radiance above 
that surface is different from the calculated theoretical value. 
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3. DATA ANALYSIS METHODS 

The investigation has utilized the radiances measured by the MSS, 
and the aerosol content measured by ground-based observations with a 
Volz photometer at the time of several ERTS overpasses. Aircraft mea- 
surements were also made of the Salton Sea/desert inherent contrast, but 
these data have not yet been received for analysis. 

3. 1 ERTS Data 

The data for the four MSS channels have been received as bulk pro- 
cessed black and white 9. 5 in positive prints and transparencies, and as 
bulk processed digital 7-track computer compatible tape, selectively ordered 
after viewing the black and white products. 

The transparencies of the first sets of data received were analyzed 
with a Jarrell-Ash microdensitometer to obtain the radiances at selected 
areas of the test sites. It was readily apparent that the photographic data 
were not accurate enough, as expected, for this investigation. On several 
transparencies, for instance, the density of the transparency over a water 
surface was greater than that of the blackest step in the calibration grey 
scale, which corresponds to zero radiance. No such obvious errors were 
found in the digital data, and all the results given in this report are based 
on the digital data. 

To extract the radiance data from the computer compatible tapes 
(CCT), a program was written to read data in prescribed geographical 
areas from the tapes on a CDC 6400. The areas of interest for analysis 
were chosen by viewing the black and white products, and selecting areas 
within the test sites free of obvious clouds, or effluents in the water. The 
location of typical selected areas, about 6 km on a side, are shown for 
each test site in Fig. 3-1. The word counts are printed out for each area, 
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(a) San Diego Test Site 


(b) Salton Sea Test Site 


Figure 3-1. Location of Analysis Areas 





and can be converted to radiance using the calibration data given in the 
ERTS Data Users Handbook. The radiance is proportional to the word 
count in each channel, with the radiance per count being given by . 0195, 

. 0157, . 0139 and . 0730 mw/cm 2 /sr for MSS 4, 5, 6 and 7 respectively. 

To analyze the Salton Sea/desert data within one of the selected 
areas, two smaller areas about 300m x 300m are chosen, one on the 
water surface and another on the desert. The water area is selected away 
from the shoreline where the radiance is essentially uniform within the 
area, and the desert area is selected for high radiance and uniformity. 

An average radiance for each area is estimated for each MSS band. 

The San Diego data are analyzed in a similar fashion, although only 
the water surfaces are of concern. The mean radiance for each MSS band 
is estimated for small areas of the ocean just off the coast, and of the 
various bays in the San Diego area. 

3. 2 Ground-Truth Measurements 

Measurements of the aerosol content were made with a Volz sun 
photometer at times of selected ERTS overpasses at the Salton Sea/desert 
site, and at every overpass at San Diego as weather permitted. 

(19) 

The Volz photometer has been widely used to measure the tur- 
bidity of the atmosphere on cloudless days. It consists of a lens, a photo- 
cell, a pivoted scale, a microammeter, and a level. When the instrument 
is level, the relative atmospheric path length (in air mass units) can be 
read directly from the pivoted scale. The instrument is then directed 
toward the sun and the radiance recorded. The solar energy enters the 
instrument through the lens whose focal length is 4. 5 cm. A diaphragm 
limits limits the angular field to about one degree. A bandpass filter 
limits the wavelength region to 0. 5 + 0. 06 micron. Before the radiation 
reaches the photocell, it is diffused by a ground glass plate mounted 
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directly behind the diaphragm. The narrow field of view used in looking 
at the small, bright source excludes nearly all scattered air light from 
the instrument, thus the aerosol optical thickness, r A , can be determined 
from 

J 

J = ~y~ exp [ -(r R + r o + r A ) m] (2-13) 

where r R and t q are the optical thicknesses due to pure Rayleigh scat- 
tering and ozone absorption, respectively, at X = 0. 5 micron and air mass, 
m, of 1. We use the values of Elterman^^ who gives t r = 0. 145 and 
t q = 0. 012. These values disagree only slightly with those used by Flowers, 
et al. who listed values as Tp = 0. 146 and t q = 0. 0092. The quantities 
J and J Q in Eq. 2-13 are the Volz photometer deflections, which are pro- 
portional to the observed solar intensity and the solar energy outside the 
atmosphere (air mass = 0), respectively. F is the sun-earth distance cor- 
rection factor. The quantity J Q is determined by measuring J over a 
wide range of air masses during periods when the aerosol optical thickness, 
r A , remains constant, and extrapolating to zero air mass. Calibrations 
of our Volz photometer several times over recent years have shown good 
repeatability. 

It was originally planned to use an Exotech radiometer to supplement 
the Volz measurements; the Exotech instrument has four channels with ap- 
proximately the same spectral response as the ERTS MSS channels. Un- 
fortunately the Exotech design does not lend itself readily to sun observa- 
tions (it was designed for downward-looking aircraft measurements), and 
no satisfactory data were obtained during the program. The first version 
of the instrument did not have baffles between the lenses and the detectors, 
so that too much scattered sky radiation was falling on the detectors. This 
problem was corrected by the manufacturer. However, it was found im- 
possible to get good repeatable readings when pointing at the sun, : without 
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significant redesign of the sighting optics and readout of the instrument, 
which was not done, and no useful data were obtained. 

3. 2. 1 Aircraft Measurements 

In order to determine the inherent contrast on the Salton Sea/ 
desert test site, two flights with the NASA C130 aircraft were planned. 
Unfortunately, the first flight scheduled for 12-12-72 was cancelled due 
to a failure of the aircraft MSS, and the next flight which took place on 
4-9-73 produced no useful data due to considerable cloud cover. Another 
flight was scheduled for, and took place on, 5-23-73 in excellent clear 
sky conditions. It was not possible, due to other priorities for the air- 
craft, to schedule another flight before the end of this program. 

The preliminary aircraft photographic data for the 5-23-73 flight 
were received and reviewed, and the digital data on computer compatible 
tapes were requested 6-26-73. Unfortunately these digital data had not 
been received at the conclusion of this program on 10-6-73, so that no 
aircraft measurements are available as ground truth. 
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4. 


RESULTS 


Significant results were obtained in analyzing the ERTS digital data 
and the ground-truth measurements. A linear relationship, as predicted 
by theory, was found to exist between the MSS radiances over water sur- 
faces and the aerosol content of the atmosphere. A linear relationship 
was also found between the contrast function (C / C 0 - 1) and the aerosol 
content. 

4. 1 Volz Data 

During the program there were twenty-two ERTS overpasses at 
each of the test sites. Six trips were made to the Salton Sea/ desert test 
site to make ground-truth measurements, and data were obtained on five 
occasions; cloud cover prevented the taking of data on just one occasion. 

It was planned to make Volz photometer measurements for each overpass 
at San Diego, but the overpasses coincided with cloud cover with surprising 
frequency, and data were obtained on only five occasions. The Salton Sea 
ground-truth measurements were always made on the shoreline at Bombay 
Beach on the east side of the sea. The San Diego observations were made 
in different locations each time, but always close to the ocean. 

The results of the measurements are given in Fig. 4-1, and show 

( 21 ) 

that in general the aerosol content was less than the value given by Elterman' 
in his 1964 model atmosphere (Elterman' s model aerosol distribution cor- 
responds to a 25 km visibility at the surface). The measurements on a few 
of the days show quite rapid changes, which are greater than the instrument 
error. They could possibly be due to the movement of undetected thin clouds 
between the sun and the observer, although all these variable measurements 
were taken in apparently completely clear sky conditions. The variations 
are probably due to the movement of air cells with different aerosol content 

I ' : 

in the lower atmosphere. 
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Fig. 4-1. Volz Data 
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4. 2 ERTS Radiance Measurements 


Although ten sets of ground -truth data were obtained at the two 

test sites, the ERTS digital data for only eight of these were received 

before the end of the program. These data, given in Tables 4-1 and 4-2, 

( 22 ) 

together with a single point provided by Fraser v ' in the Atlantic (see 
Table 4-3), are used to investigate the radiance-aerosol content relation- 
ship, and the contrast -aerosol content relationship. The radiance values 
for MSS 7 are not listed in the tables, since they are low and inaccurate 
over water surfaces, being close to the instrument noise level (see Section 
4. 3. 2). In addition, the radiance in this channel is influenced by the 
variable atmospheric water vapor content. 

4. 2. 1 Radiance-Aerosol Relationship (Water Surface) 

The investigation of the radiance -aerosol relationship over water 
surfaces is based on eight data points, four over the Pacific at San Diego, 
three over the Salton Sea, and one over the Atlantic off the coast of North 
Africa, where measurements of the aerosol content are available. The 
radiance values are normalized to a sun angle of y = 0. 45, using a cor- 
rection based on the theoretical variation with sun angle given in Fig. 2-3; 
it is assumed that the correction is the same for each channel. The cor- 
rection is greater than 5% only for the three points with high sun angle. 

The normalized radiances for all four MSS channels are plotted 
against the aerosol content in Fig. 4-2. It is seen that, as predicted by 
theory, a linear relationship exists for all four channels, being best for 
MSS 5 and MSS 6. MSS 6 shows excellent agreement with the theoretical 
relationship shown in Fig. 2-3, suggesting that the model is very reason- 
able, particularly for zero aerosol content (i. e. , a pure molecular atmo- 
sphere). The results for MSS 7 are given only to illustrate the low radiance 
values discussed above, and even though a linear relationship is indicated, 
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TABLE 4-1. 


Radiance (mw/cm YMin/sr) Data for San Diego Test Site 


Date 

Sun 

Angle 

(#*) 

Volz 

Aerosol 

Content 

MSS 

Channel 

La Jolla 
Ocean 
Radiance 

Mission Beach 
Ocean 
Radiance 

Point Loma 
Ocean 
Radiance 

Mission 

Bay 

Radiance 

San Diego 
Bay 

Radiance 

11-25-72 

0. 52 

0.42N 

4 

2. 92 

3.32 

3. 51 

3.32 

3. 51 



(Pt. Loma) 

5 

0.95 

1. 10 

1. 10 

1. 10 

1.26 




6 

0.42 

0. 47 

0. 44 

0. 42 

0. 51 

12-13-72 

0.47 

1. 01 N 

4 

3.32 

Cloud 

Cloud 

Cloud 

Cloud 



(La Jolla) 

5 

1.42 








6 

0. 69 





12-31-72 

0.45 

0. 39 N 

4 

3. 12 

3. 12 

3. 12 

3.22 

3. 12 



(La Jolla) 

5 

0.95 

1.03 

0.95 

1. 18 

1. 18 




6 

0.42 

0.42 

0.42 

0. 55 

0. 55 

1-18-72 

0. 47 

0. 51 N 

4 

2.92 

3. 42 

3. 51 

3. 61 

3. 51 



(La, Jolla) 

5 

1. 00 

1.26 

1.26 

1.42 

1.42 




6 

0.42 

0. 55 

0. 49 

0. 62 

0. 55 


2 

TABLE 4-2. Radiance (mw/cm /um/sr) Data for Salton Sea/Desert Test Site 


Date 

Sun 

Angle 

GO 

Volz 

Aerosol 

Content 

Inferred 

Aerosol 

Content 

MSS 

Channel 

Radiance 

Salton 

Sea 

Desert 

Apparent 

Contrast 

C R 

12-12-72 

0.47 

0. 54 N 

0. 64 N 

4 

2.90 

9.20 

2.2 





5 

1. 10 

8.30 

6. 5 





6 

0. 55 

6. 50 

10.7 

4-17-73 

0.82 

1. 19 N 

1. 19 N 

4 

4. 68 

15.8 

2.4 





5 

2.20 

15.9 

6.2 





6 

1. 11 

12.5 

10.3 

5-5-73 

0.87 

0. 77 N 

1. 46 N 

4 

5.23 

17.5 

2. 1 


(Effluents present in water) 


5 

2. 68 

16.6 

4. 1 





6 

1. 53 

13.4 

6.4 

5-23-73 

0.88 

1.11N 

1.05N 

4 

4.87 

18. 1 

2.7 





5 

2. 42 

18.3 

6.6 





6 

1. 21 

14.3 

10.9 

8-26-72 

0.82 


0.98N 

4 

4.20 

14.6 

2.5 





5 

2.05 

14. 5 

6. 1 





6 

1.04 

11. 1 

9.7 

9-13-72 

0. 77 


1. 18 N 

4 

3.70 

14.0 

2.8 





5 

1. 80 

13.9 

6.7 





6 

1.07 

11. 1 

9.4 



Table 4-2 continued. . 


Date 

Sun 

Angle 

in) 

Volz 

Aerosol 

Content 

Inferred 

Aerosol 

Content 

MSS 

Channel 

Salton 

Sea 

Radiance 

Desert 

Apparent 

Contrast 

C R 

10-1-72 

0. 72 


1. 12 N 

4 

3.90 

14.8 

2.8 





5 

1.70 

13.4 

6.7 





6 

0. 97 

10. 3 

9.6 

11-6-72 

0. 57 


0. 91 N 

4 

3.90 

12.9 

2.3 





5 

2.00 

11.8 

4.9 





6 

1.20 

9. 1 

6.6 

11-24-72 

0. 52 


0.73N 

4 

3. 12 

10. 3 

2.3 





5 

1.29 

9.8 

6.6 





6 

0. 60 

8.0 

12. 5 

12-30-72 

0.45 


0. 52 N 

4 

2. 83 

9.45 

2.4 





5 

1. 10 

8.90 

7.1 





6 

0.49 

7. 10 

13.6 

1-17-73 

0.47 


0. 53 N 

4 

3. 12 

9.35 

2.0 





5 

1. 11 

8.98 

6. i 





6 

0. 55 

6.92 

11.5 

3-30-73 

0.75 


1. 08 N 

4 

4. 38 

15.2 

2.5 





5 

2. 12 

15. 1 

6. 1 





6 

1. 01 

11. 9 

12.2 



TABLE 4-3. 


2 

Radiance (raw/cm //jm/sr) Data for 
Atlantic Ocean (21°N, 17°W) (22) 


Date 

Sun 

Angle 

00 

Volz 

Aerosol 

Content 

MSS 

Channel 

Radiance 

8-9-72 

0. 86 

2.33 N 

4 

5. 52 




5 

3. 18 




6 

2.04 


30 


Radiance (mw/cm V^m/sr) 



Fig. 4-2. Radiance vs. Aerosol Content Over Water Surfaces 
(Normalized to sun angle = . 45) 
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this channel is not considered useful for determining the aerosol content. 
The results for MSS 4 show more scatter of points about the mean curve 
than for the other channels. The reason for the larger scatter in this 
channel is not readily obvious, since the radiance errors are not signifi- 
cantly different from the other channels (see Section 4. 3. 2). Possibly 
this channel is more sensitive to changes in the water reflectance due to 
the presence of effluents, chlorophyll, or other suspended material in the 
water, although on the one occasion when effluents were obviously present 
in the Salton Sea, the radiance in MSS 4 was closer to the mean curve than 
in the other channels. 

Effluents were clearly apparent on 5-5-73 when ground -truth mea- 
surements were made, and on subsequent imagery of the Salton Sea during 
the summer months. The effluents, which are presumably due to irriga- 
tion run-off, were not observed during the winter months, when there is 
more rain and it is cooler, thus requiring less irrigation. These effluents 
result in variable high radiances over the water surface. To analyze the 
data for 5-5-73, an area exhibiting the lowest values of radiance was 
selected, but as seen in Fig. 4-2, the radiances appear too high for the 
measured aerosol content, especially for MSS 5 and 6. The effluent pattern 
is more apparent in the photographic images at shorter wavelengths. How- 
ever, the data point for MSS 4 shows less deviation from the previous data 
than at the longer wavelengths. In addition, the spectral variation of ra- 
diance has the same basic shape as other data. Thus, the spectral be- 
havior of this effluent, and its effect on aerosol observations is not com- 
pletely consistent, and further studies of this type of data are required. 
However, it is clear that care should be taken in using data from bodies 
of water where effluents occur on an intermittent basis. The data may be 
readily screened by visual examination of the photographic images. 
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4.2.2 Radiance-Aerosol Relationship (Desert Surface) 

The theory discussed in Section 2 showed that the radiance over 
high albedo surfaces is not very sensitive to atmospheric aerosol content 
changes. The data at the Salton Sea/desert test site allows this conclusion 
to be investigated over the desert, which has a high albedo (~ 0. 3). 

Since ground truth are available for only four ERTS overpasses, the 
data for several other ERTS overpasses are utilized with help of the ra- 
diance-aerosol relationship in Fig. 4-2. This linear relationship is assumed 
to be correct, and is used to determine the aerosol content for twelve over- 
passes from the Salton Sea radiance for MSS 6. The desert radiances, given 
in Table 4-2, are plotted in Fig. 4-3, against these inferred aerosol contents, 
also listed in Table 4-2. The data, uncorrected for sun angle, appear to 
show a good linear relationship. However, after normalizing the radiances 
to a sun angle of .45, using the theoretical variation for a Lambertian 
albedo of 0. 3 (actually 0. 4 or 0. 2 which show approximately the same var- 
iation) given by Plass and Kattawar, v ' the radiance clearly shows no obvious 
dependence on the aerosol content, as predicted by theory. It is realized 
that the desert is not a true Lambertian surface, with the backscattering 
being greater than the forward scattering. ' ' However, since it is close 

to being a diffuse reflector, and no significant variability of albedo with 
sun angle has been found, it is believed that the above conclusions are 
not affected. 

4. 2. 3 Analysis of Potential Problem Areas 

Two potential problem areas for this program were identified in 
Section 2. 3. The first one, sun glitter, may be discussed on the basis of 
the preceding radiance-aerosol content analysis, and the second, surface 
reflectance gradients, should be studied before analyzing the contrast - 
aerosol content relationship. 
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Fig. 4-3. Desert Radiance vs. Aerosol Content for MSS 6 


34 


4. 2. 3. 1 Sun Glitter 


No evidence of sun glitter has been observed in the ERTS data 
received during the program. (Sun glitter was observed on the photo- 
graphic data from the NASA aircraft overflight at 2000 ft. altitude on 
5-23-73, but was not identifiable on the ERTS data for the same 
time). All of the data at the San Diego and Salton Sea/desert test sites 
were obtained over relatively calm seas, so that the sun glitter might be 
expected to be minimal. However, the one data point in Fig. 4-2 for the 
Atlantic Ocean was obtained over rough seas, with seven foot waves and 
winds of about 20 knots being reported. Presumably, sun glitter should 
be apparent under these conditions. However, it is seen in Fig. 4-2 that 
this point does not significantly deviate from the linear relationships for 
the San Diego and Salton Sea data. Thus it is believed that nadir observa- 
tions will not be significantly affected by sun glitter, although more rough 
sea data are needed to draw any firm conclusion. 

It was suggested in Section 2. 3. 1 that observations at two wave- 
lengths might be used to eliminate a sun glitter problem if it should occur. 

In order to investigate this possibility the spectral variations of the upwelling 
radiance over water surfaces have been plotted. Figure 4-4 shows the 
results for several ERTS overpasses at the Salton Sea. These values are 
not normalized, and higher radiances are expected for higher sun elevations 
for a given aerosol content. Within the errors of the MSS radiance mea- 
surements, there are no obvious spectral differences from one overpass 
to the next, at least at the MSS wavelengths. However, there appears to 
be a tendency at sun elevations below about 50° for the radiance decrease 
from MSS 5 to MSS 6 to be less than at high sun elevations. 



Radiance (mw/cmV^Wsr) 



Wavelength (pm) 


Fig. 4-4. Salton Sea Water Radiance vs. Wavelength 






The data for the ocean and bays at the San Diego test site are given 
in Fig. 4-5, and show the same basic variations as the Salton Sea mea- 
surements. Also given in Fig. 4-5 are the data for the Atlantic Ocean 
provided by Fraser. These radiances, which are averaged over a 
20 x 20 km area show less decrease with wavelength than the San Diego 
and Salton Sea data. The Atlantic Ocean was rough, whereas the San Diego 
data, and probably all of the Salton Sea data, were for relatively calm 
surfaces, without whitecaps. If the difference in the Atlantic data is due 
to the surface roughness, then it is unlikely that a two wavelength method 
could be used to eliminate a sun glitter problem. However, with only one 
set of rough sea data, and the errors in the measurements, no conclusion 
can be reached at the present time. 

4. 2. 3. 2 Surface Reflectance Gradients 

Examination of the digital data in the San Diego coastal regions and 
at the Salton Sea shows that the radiance over the water surface is higher 
immediately at the water -land boundary, but rapidly decreases and reaches 
a uniform value within about 1300 ft of the boundary. Thus, the effects of 
scattering from the adjacent high albedo land appears to be negligible beyond 
about 1300 ft from the shoreline. The water radiance data used for the 
preceding radiance-aerosol study and the following contrast -aerosol analysis 
were all obtained beyond 1300 ft from the shoreline, where the radiance was 
uniform. 

4. 2. 4 Contrast -Aerosol Relationship 

The Salton Sea/desert test site is used to investigate the contrast - 
aerosol content relationship discussed theoretically in Section 2. 2 To de- 
termine the function f(r) given in Eq. (2-12), and plotted in Figs. 2-4, 2-5, 
and 2-6, the value of A', the water surface reflectivity, and C Q , the 
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Fig. 4 -5(a). San Diego and Atlantic Ocean 

Water Radiance vs. Wavelength 


IT' 
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inherent contrast, must be known as a function of wavelength and sun 
angle; C R , the apparent contrast, is calculated directly from the ERTS 
digital data, and is listed in Table 4-2. Since the NASA aircraft MSS data 
are not available to determine C Q , and A' is not known, certain ap- 
proximations must be made. A' is assumed to be independent of sun 
angle and wavelength so that the function (C q /C r - 1) (= f(r)/A') is 
examined. C is estimated by assuming a reflectivity of 0. 30 for the 
desert and 0. 02 for the Salton Sea (based on albedos measured by Griggs' ') 
giving a value of 14. 0 for C Q . It is further assumed that C Q is inde- 
pendent of wavelength and sun angle. 

Ground -truth measurements are available for only four ERTS 
overpasses at the Salton Sea/desert test site, and on one of those occasions 
effluents were present in the water, so that only three data points can be 
used. However, the data from several other ERTS overpasses are uti- 
lized with the help of the radiance-aerosol content relationship in Fig. 4-2. 
This linear relationship is assumed to be correct, and is used to deter- 
mine the aerosol content from the Salton Sea radiance for MSS 6. These 
values are then treated as ground-truth values of the aerosol content. 

The values of (C /C_, - 1) are calculated for each overpass, and 
normalized to a sun angle of y= 0.45, using the theoretical relationships 
in Figs. 2-4, 2-5, and 2-6. The variation of (C q /C r - 1) with sun angle 
is due to the variation of C R alone. The function (C Q / C R - 1) is plotted 
in Fig. 4-6 against the aerosol content, as defined in Section 4. 2. 1, for 
MSS 4, 5 and 6. The intercept of the abscissa is determined for the hypo- 
thetical case of no atmosphere when C Q = C R , and is given by the nega- 
tive value of (t r /t^)N where r R is the Rayleigh optical thickness, 
is the Elterman 1964 model aerosol optical thickness, both for the center 
wavelength of the particular MSS channel, and N is the Elterman 1964 
model aerosol content. ; : 
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Fig. 4-6. (C o /Cp - 1) vs. Aerosol Content 
(Normalized to sun angle U = . 45) 
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Figure 4-6 shows that good linear relationships exist for MSS 4, 

5 and 6; MSS 7 is not considered due to its low inaccurate radiance values 
over water. The results show that the contrast function becomes more 
sensitive to aerosol changes as the wavelength decreases. This is due to 
the increased atmospheric optical thickness at the shorter wavelengths. 

It should be noted that, since the radiance over the desert surface is not 
very sensitive to aerosol changes (see Section 4. 2. 2), the contrast function 
depends mainly on the radiance over the water surface. 

When the instrument errors, discussed in Section 4. 3, are con- 
sidered, it is seen that MSS 6 is not sensitive enough, but that MSS 4 and 
5 could be used for this technique and that they have comparable accuracy. 
If the mean straight lines are correct, then the instrument error in both 
MSS 4 and 5 would produce about + 10% errors for the normal aerosol 
content. The scatter of data points about the mean line in Fig. 4-5 is due 
to measurement errors in the desert radiance, and due to errors in the 
assumptions concerning the contrast dependence on wavelength and sun 
angle, which can be examined only with independent measurements of these 
properties. 


4. 3 Error Analysis 

The results of the analyses of the radiance -aerosol content and 
contrast -aerosol content relationships contain uncertainties due to inherent 
instrument errors, and due to the assumptions necessary in the data inter- 
pretation. These uncertainties are discussed below. 

4. 3. 1 Volz Photometer Errors 

Flowers et al, in comparing the Volz photometer with a stand- 
ard photometer, found that values of aerosol optical thickness greater than 
about . 240 (i. e. , an aerosol content of 1. 13 N), can be measured with an 
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accuracy of i5%, and for an optical thickness near . 100 (.47 N), an accuracy 
of about ± 10% applies. It is presumed that similar errors apply to the 
measurements made in this program with our Volz photometer, and are 
indicated in Figs. 4-1 and 4-2. 

A further error is possible in using the Volz data in this investi- 
gation due to the fact that the Volz measurement was never made directly 
over the water surface where the ERTS radiance measurement is made. 

It has been assumed that the atmosphere has been homogeneous over the 
area including the Volz and the radiance measurements. This may not 
always be a good assumption since the Volz data in Fig. 4-1, on occasion, 
show variations over a short time period, greater than the instrument 
error, suggesting the movement of air cells with different aerosol con- 
tent. This type of error cannot be readily estimated, but could be 
checked in a future program, using several photometers in the test site 
area. 

4. 3. 2 ERTS Radiance Errors 

Fraser has made calculations of the noise equivalent radiance 
for each of the MSS channels based on preflight calibrations. No inflight 
calibrations have been reported, and it is assumed that the preflight ones 
have not changed during this program. Fraser's results are reproduced 
in Table 4-4, with the values of the full scale radiance, and NER in ra- 
diance units added for clarity. If the radiance-aerosol content relationship 
in Fig. 4-2 is assumed to be correct, and the reduction in noise achieved by 
averaging over a 300m x 300m area (approximately 25 resolution elements) 
is taken into account, then the NER results in about a ± 10% error in deter- 
mining the aerosol content in the useful channels MSS 4, 5 and 6. 

If the radiance -aerosol content relationship in Fig. 4-2 is assumed 
to be correct, then the NER results in a ± 10% error in determining the 
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aerosol content in the useful channels MSS 4, 5 and 6. A similar uncer- 
tainty of i 10% in the aerosol content also results from the NER in the 
contrast -aerosol content relationships for MSS 4 and 5 shown in Fig. 4-6. 

With the limited set of data obtained during this program, certain 
assumptions were necessary regarding the reflectance of water and the 
desert as a function of wavelength and sun angle. These assumptions, 
as discussed in Section 4. 2, were based on available measurements and 
theoretical models, and the uncertainties introduced cannot be readily 
estimated. The relatively small scatter of data in the linear relation- 
ships suggest that the assumptions are not unreasonable. 


TABLE 4-4 

The noise equivalent radiance (NER) for a single resolution element 
of MSS for ocean observations. Also the ratio of the standard deviation 
in sensor response to atmospheric turbidity changes to the NER. The 
mean sensor output (V), the standard deviation (a), and NER are 
expressed in percentage of full scale sensor response (after Fraser ). 


Spectral 
Band 
in Mm 

Surface 

Reflectivity 

V 

a 

SNR 

V/SNR 

c/NER 

Full Scale 
Radiance 
(mw/cmVsr) 

NER 2 
mw/ cm/ 
ym/sr 

0. 5-0. 6 

0. 02 

12 

1.4 

22 

0.6 

2.5 

2.48 

.15 

(MSS 4) 

0. 10 

23 

1.1 

40 

0.6 

1.9 



0. 6-0. 7 

0.02 

8 

1.4 

9 

0.9 

1.5 

2.00 

.18 

(MSS 5) 

0.10 

22 

1.0 

22 

1.0 

1.0 



0. 7-0. 8 

0.02 

7 

1.0 

11 

0.6 

1.6 

1.76 

. 11 

(MSS 6) 

0. 10 

22 

1.0 

25 

0.9 

1.1 



0.8-1. 1 

0. 02 

4 

0.5 

3 

1.0 

0.5 

4. 60 

. 15 

(MSS 7) 

0.10 . 

10 

0.4 

10 

0.9 

0.4 
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5. 


CONCLUSIONS 


Significant results, relating the radiance over water surfaces to 
the atmospheric aerosol content, have been obtained. The results indicate 
that the MSS channels 4, 5 and 6 centered at 0. 55, 0. 65 and 0. 75 fim 
have comparable sensitivity, and that the aerosol content can be deter- 
mined within + 10% with the assumed measurement errors of the MSS. 

The fourth channel, MSS 7, is not useful for aerosol determination due 
to the water radiance values for this channel generally being less than 
the instrument noise. The accuracy of the aerosol content measurement 

could be increased by using an instrument specifically designed for this 

( 27 ) 

purpose. In an independent study v ' we designed a simple instrument in 
which the radiance could be measured to ± 1% accuracy resulting in a 
+ 1.5% error in the aerosol content. 

This radiance -aerosol content relationship can provide a basis for 
monitoring the atmospheric aerosol content on a global basis, allowing a 
base-line value of the global burden of aerosols to be established. This 
base-line could be established more rapidly from satellite measurements 
than from a network of ground-based observations, and probably with 
considerable cost savings. In addition, this technique could provide a 
method for monitoring the particulate emissions of the SST's, which are 
of concern to the Department of Transportation, by making observations 
in the vicinity of flight corridors, such as over the North Atlantic. It 
should be possible to look at the ocean through the flight corridor and 
alongside it, to measure the difference due to the SST’s in the aerosol 
content. 

Further studies of the radiance -aerosol content relationship should 
be made to establish the global applicability of the results, and to confirm 
that the effects of sun glitter are minimal as indicated in this program. 
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These studies could utilize future ERTS-1 and ERTS-B data at the existing 
test sites. Cooperation with other agencies such as EPA and NOAA, who 
make Volz observations, would allow the studies to be expanded to national 
and global scales. ERTS-1 data obtained over the USA since August 1972 
could also be analyzed in conjunction with aerosol data from the turbidity 
network of Volz photometers operated by EPA. 

The contrast -aerosol content investigation showed useful linear 
relationships in MSS channels 4 and 5, allowing the aerosol content to be 
determined within ± 10%. MSS 7 is not useful due to the low accuracy in 
the water radiance, and MSS 6 is found to be too insensitive. These results 
rely on several assumptions due to the lack of ground-truth data, but do 
serve to indicate which channels are most useful. 

Future studies of this contrast relationship should use a different 
target where the radiances over both surfaces are affected significantly 
by aerosol changes. Possibly the Salton Sea/Imperial Valley area, or the 
San Diego bay/city area could be used, but several simultaneous Volz ob- 
servations would be needed to assure homogeneity of the aerosol content 
over the target. The data from the EPA turbidity network should be exa- 
mined for suitable target areas to investigate this relationship with existing 
ERTS data. 
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